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Abstract — Image resizing is necessary when you need to 
increase or decrease the total number of pixels, whereas 
remapping can occur under a wide variety of scenarios. If the 
same image resizing or remapping is performed, the results can 
vary significantly depending on the interpolation algorithm. In 
this paper, image magnification performance of different 
interpolation algorithms was compared from subjective and 
objective aspects. The experimental results give the guidance 
for the user to choose a suitable algorithm to achieve optimum 
results according to different applications. 

Index Terms — Image remapping, Image resizing, 
Interpolation algorithm. 


I. Introduction 

Image resizing is one of the basic image operations, and 
is widely used in many applications [1]. The image 
magnification is a conversion process from a low resolution 
image to a high resolution image. The image magnification 
is essentially image interpolation process. There have been a 
lot of practical image magnification methods which have 
their own characteristics, advantages and disadvantages. The 
choice of the different image magnification algorithm 
directly affects enlarged image quality, so it is important to 
find a suitable algorithm to improve the image magnification 
quality. In current, the commonly used image magnification 
method is interpolation method. In the paper image 
magnification performance of different interpolation 
algorithms was compared from subjective and objective 
aspects to give the guidance for the user to choose a suitable 
algorithm to achieve optimum results according to different 
application. 


enlarged, for example if an image is zoomed to 400% shown 
in Fig. 1., the color values of original 4 adjacent pixels 
marked A, B, C, and D in (a) were filled in the new A, B, C, 
and D locations in (b) accordance with the magnification 
factor. But there are a large number of pixels which values 
are unknown between A, B, C, and D, such as P. So the 
values of these pixels should be calculated through 
interpolating estimation. 
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Figure 1 Schematic diagram of Image Magnification 

III. COMMON INTERPOLATION ALGORITHM 


A. Nearest Neighbor Interpolation 

In nearest neighbor interpolation algorithm, the position 
of pixel P in the magnified image is converted into the 
original image, and the distance between P and its neighbor 
points A, B, C and D were calculated. Then the color values 
of pixel P was set as the values of the pixel which was nearest 
toP [2]. 


II. IMAGE AMPLIFICATION PRINCIPLE 

The basic principle of image magnification is to increase 
the image pixel number, so a low resolution image is 
converted to a high resolution image. When a small image is 
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In Figure. 2, suppose (i, j ) , (i, j + 1) , (i+ 1, j ) and ( 
i+l,j+ 1) are the 4-neighbor points, and there values are f( i, 
j) , f( i, j + 1) , f( i+ 1, j ) and f( i+ l,j + 1). The distance 
between (u,v) and (i, j), (i, j + 1), (i+ 1, j) and (i+ l,j + 1) 
were calculated, then the values of (u,v) was set as the value 
of the point which is nearest to (u,v). [3] 
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Figure 2 Diagram of Nearest Neighbor Interpolation Algorithm 

A. Bilinear Interpolation 

In bilinear interpolation algorithm [4], the position of 
pixel P in the magnified image was converted into the 
original image, then the Influence of the four pixel points 
A,B, C and D was calculate. The nearer distance to the point 
P, the value is greater, which indicates the greater effects. 
The diagram of bilinear interpolation is shown in Figure. 3 
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Figure 3 Diagram of Bilinear Interpolation Algorithm 


Suppose the coordinate of A, B, C and D are (i, j), (i, j + 
1), (i+ 1, j) and (i+ 1, j + 1). The coordinate of P is (u,v). 
The bilinear interpolation algorithm consists of three steps: 

Step 1: To calculate the influence of A and B and denote it as 

E. f (i, j+v) = [f (i, j+l)-f (i, j)] v+f (i, j) 

Step 2: To calculate the influence of C and D and denote itas 

F. f (i +1, j+v) = f (i +1, j+1 )-f (i+1, j)] v+f (i+1, j) 

Step 3: To calculate the influence of E and F and denote it as 
P.f(i+u,j+v)=( 1 -u)( 1 - v)f(i,j)-( 1 -u)vf(i,j+1 )+u( 1 - v)f(i+1 ,j)+uv 

f(i+lj+l) 

B. Bicubic Interpolation 

Bicubic interpolation is similar to bilinear interpolation 
algorithm. For the unknown pixel P in amplified image, its 
influence sphere is expanded to its 16 adjacent pixels, then 
the color value of P is calculated by these 16 pixels according 
to their distance to P. The diagram of bicubic interpolation 
algorithm is shown in Figure. 4. 



Figure 4 Diagram of Bicubic Interpolation Algorithm 


Compared with bilinear interpolation algorithm, 
bicubic interpolation algorithm extends the influence with 


more points, and uses advanced interpolation algorithm. In 
the horizontal direction, to calculate the value of e between 

two points A and B, it is need to use four pixels 
value of A, B , A-l, B+l and get a smooth curve through a 
nonlinear calculation. The diagram of nonlinear 
interpolation algorithm is shown in Figure 5. 
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Figure 5 Diagram of nonlinear interpolation algorithm 

On the basis of calculating the lateral interpolation 
effects, the operation as mentioned above is extended into 
two dimensional space to calculate the longitudinal 
interpolation effects. 

C. Cubic B-Spline 

In order to get better image quality, many new algorithm 
shave been proposed, they use more complex improved 
interpolation method such as B spline [5] and Mitchell 
interpolation algorithm. Their objective is to make the 
interpolation curve more smoothly and the image edge more 
perfect. The function of cubic B-spline is as follow [2]. 

yrx)=2fr=-*5fr.TH- iCOTO*) w 

Three order B-spline function is as follows: 
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IV. EXPERIMENT RESULT 

In order to test the performance of the non-linear 
interpolation algorithms, We have developed a MATLAB 
2009 [6] implementation on Intel(R) Pentium(R) Dual E2160 
2.00GHz CPU, 2.0GB RAM computer. To compare the 
performance of the different algorithm, a test image was 


< X < X i+ 1 

,x i+l <x<x i+2 

(5) 

,x i+ 2 <x<x i+3 

everywhere 
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shrunken to half of the original size. And then the shrunken 
image was expanded to its original size by using different 
algorithms. At last the enlarged image quality was compared 
from two aspects: Subjective and Objective evaluation. 


Figure 6 Computational Time of different interpolation algorithms 
B. Subjective Evaluation 


A. Objective Evaluation 


In this paper, the peak signal-to-noise ratio (PSNR) was 
used to compare the reconstructed image with the original 
image in image compression. 


PSNR = 20 log :o 


'MAX r 

— 

.vmseJ 



In the experiment, we selected LENA Image; image to test, 
the results is shown in Table I. 


Table I PSNR for Different Interpolation Algorithms. 


Interpolation 

Type 

PSNR(dB) LENA 
Image 

Nearest Neighbour 

26.05 

Bilinear 

27.183 

Bicubic 

27.184 

B-Spline 

27.150 


PSNR can reflect the magnified image quality; the 
larger PSNR means the higher image quality. It can be seen 
from the above table that the SNR between the original image 
and the magnified image get from the nearest neighbour 
interpolation is the minimum, and bilinear interpolation 
algorithm is much larger, bicubic interpolation is the largest, 
and the B-Spline larger than the nearest neighbor method but 
less than the bilinear interpolation algorithm and bicubic 
interpolation method. The complexity of different 
interpolation algorithm is different, and their computational 
time also has bigger difference. The computational time is 
shown in Figure 6. 


Computational Time 



Image is used to be watch by people, so it is more suitable 
evaluate the magnified image quality with subjective 
methods. In this paper, 8 observers were invited to observe 
magnified images and give evaluation. At last the image 
scoring average and used to evaluate the subjective image 
quality. The subjective evaluations with different 
interpolation methods are shown in Table II. 


Table II Subjective Evaluation with Different Interpolation 

Methods 


Interpolation 

Type 

Subjective 

Feelings 

Image 

Contour 

Overall 

Evaluation 

Processing 

Time(secs) 

Nearest Neighbour 

Mosaic 

phenomenon 

Not Clear 

Worst 

5 

Bilinear 

Blur. not 
sharp 

Not Clear. Senate 
phenomenon 

Poor 

6 

Bicubic 

Sharper 

Serrate phenomenon 
improved 

Better 

8 

B-Spline 

Relatively 
clear and 
sharp 

Edge becomes clear, 
Serrate phenomenon 
disappeared 

Good 

17 




a)Original Image b) Contract Image 



c) Nearest Neighbour d) Bilinear Interpolation 

Interpolation 



e) Bicubic Interpolation 


f) B-Spline Interpolation 


Figure 7 Different non-linear interpolation algorithms result 
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V. Conclusion 

Nearest Neighbour Interpolation algorithm is the most 
simple and fast algorithm. It has the advantages of fast speed, 
but it can bring significant distortion and it will appear 
mosaic and saw tooth phenomenon. 
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Bilinear Interpolation method is more complex than the 
nearest neighbour method, and so it has larger calculation. It 
has no gray discontinuity defects and has satisfactory results. 
This method has a low pass filtering properties, so that the 
high frequency component is faded and the image contour 
has some degree of fuzzy. As the results of bilinear 
interpolation algorithm is continuous, so the visual effect is 
better than the nearest point interpolation algorithm, but the 
operation speed is slightly slower. Because of this, if someone 
pays more attention to speed and good visual effect, this is a 
good compromise. 

Bicubic Interpolation algorithm can get relatively clear 
picture quality, but it needs larger amount of calculation. 
This algorithm now is most commonly used in many image 
processing software such as Photoshop, After Effects, Avid 
and Final Cut Pro etc. 

The magnified image by using of Cubic B-spline 
Interpolation algorithm is smooth and has no obvious saw 
tooth phenomenon. But in the interpolation process it takes 
the form of a low pass filter and suppress high frequency 
component in varying degrees. When the amplification 
factor is higher, it will cause the edge fuzzy and false 
artificial traces. 

The experimental results give the guidance for the user 
to choose a suitable algorithm to achieve optimum results 
according to different application. 
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